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)—{ Abstract 



The effect of an inhomogeneous magnetic field which varies inversely as distance on the ground 
state energy level of graphene is studied. In this work, we analytically show that graphene under 
the influence of a magnetic field arising from a straight long current-carrying wire ( proportional to 
i-^ the magnetic field from carbon nanotubes and nanowires) exhibits zero energy solutions. We find 

^H that contrary to the case of a uniform magnetic field for which the zero energy modes show the 

O localization of electrons entirely on just one sublattice corresponding to single valley Hamiltonian, 

zero energy solutions in this case reveal that the probability for the electrons to be on the both 
sublattices, say A and B, are the same. 
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^ 1 Introduction 

> 

• ^H Graphene, a single layer of graphite, was isolated for the first time in 2004 [1]. The carbons atoms 

/\ in graphene are arranged into a honeycomb structure which is consistent of the two inequivalent 

^ triangular sublattices, say A and B [2]. Electrons in graphene can hop to the nearest neighbours 

atoms which leads to the formation of the two energy bands, each containing the same number of 
states [3] and touching each other at the two inequivalent points called Dirac points, say K^ and K~ . 
Around these points the energy dispersion relation of graphene is linear in momentum which implies 
that its low energy excitations mimic the ultra relativistic massless particles. Thus, the low energy 
excitations of graphene are described by the following Dirac-like equation: 

H = vpo-.Tp (1) 

where vp is the Fermi velocity and cr={ax,cry) is the Pauli matrices vector with at, i = x,y,z, the i 
Pauli matrix. The above equation implies that the electrons in graphene behave as massless charged 



'Electronic mail: rezakord@ipm.com 
'Electronic mail: DariuslillO@gmail.com 



Dirac ferniions confined in a 2D space, an interesting feature that real particles do not exhibit because 
all the massless elementary particles happen to be electrically neutral. These niassless electrons shows 
peculiar properties which massive relativistic carriers do not exhibit [4] . In fact the first experimental 
evidence that revealed the charged carriers in graphene mimic massless electrons was an unusual 
quantum Hall effect reported in 2005 [5] . In spite of this fact that charge carriers in graphene exhibit 
a four fold degeneracy (which comes from the real spin of electrons in addition to another factor of 
two, due to the equal contributions of the K- valleys, i.e. K"*" and K^) we see that in experiments 
instead of quantization of the Hall conductivity in multiples of 

(2) 



axy = 4n— ne{...,-2, 


-1,0,1,2,...} 


uctivity is: 




1 e2 e^ 
-.. = 4(n + -)-=4.-, 


u = {n+l] 



(3) 

which shows that the integer quantum Hall effect (IQHE), appears in half- integers. It is also ob- 
served that unlike to the quantum Hall effect for 2D conventional systems which appears in the strong 
magnetic field and low temperature limit, the IQHE in graphene can be observed even at the room 
temperature [6]. This is because of ultra-relativistic nature of its charge carriers which mimic the 
massless Dirac fermions. These massless charge carriers as we'll show later, contrary to the conven- 
tional 2D systems, show interesting results under the influence of a uniform magnetic field. We, before 
discussing the effect of a constant magnetic field B = i?^ perpendicular to the graphene's plane, note 
that even for conventional 2D systems, the periodic potential due to the host lattice is of no relevance 
to the quantum Hall problem because the size of the electron wave packet in a magnetic field is much 
larger than the lattice period. The periodic potential due to the lattice is, therefore, neglected in 
studies of the quantum Hall effect, however if one considers to calculate the Landau levels based on 
the tight-binding model, the commensurability problem between the magnetic flux and lattice unit 
cell is needed to be considered. This problem is known to inevitably occur in the two dimensional 
electron system [7-9]. However, interestingly for graphene the periodic potential of the honeycomb 
lattice is already built-in and therefore it is counted in the massless Dirac-like equation. Thus, we do 
not really need to incorporate explicitly the periodic potential term into the Dirac equation. 
Now, in order to obtain the energy spectrum of graphene in the presence of a uniform magnetic field 
which is considered to be perpendicular to the garphene's plane, by choosing the symmetric gauge 
and taking the units such that c = 1, the single valley Hamiltonian of graphene can be written as: 

„ / Ux-iUy\ 

^ = "Hn. + m, J' (^) 

where 

n = -zM9.-i^,9, + ^), /^ = ^ (5) 

Then, one may write the equation (4) in the form of the following eigenvalue equation: 

{ihdt + ifiVpCTxIlx + i^VFcryIly)'^j^+{r, t) = 0. (6) 

Next, multiplying the above equation by o"^ (the z-component of Pauli matrix) gives: 

{ihj^dt + invFj^Hx + ihvF7^Uy)^j^+ (r, t) = 0, (7) 



with the 7 matrices as: 
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-i 
-i 
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Here in order to solve the equation (7), we spht the 2-spinor ^k:+ ^'^^'^ ^^^ sublattice parts: 



(8) 



^K+(r,t) 






(9) 



which inserting it into the equation (7) leads us to the following expression: 



•f 

-^ 



S 
D 






where S and D are given by: 



2/2' 



S = -vpidx - WP?) + ivpidy + ^), 



D = vpidc - ^) + ivpidy + ^ 



(10) 



(11) 



(12) 



There is no need to say that from equation (10) one can obtain the following second order differential 
equations: 



If = —SDip, 
i) = -DSip, 

which by introducing the following dimensionless quantities: 

X , ~ y 

X -^X= J, y -^y= J, 



(13) 
(14) 



(15) 



we can solve the them with respect to cp and ip. In order to do so, we first need to make the following 
ansatz: 



ip{r,t) = e-'^^/^ip{r), 



and 



V^(r,t)=e-^*/^^(r), 
which plugging them in the equations (13) and (14) yields: 



2„,2 



E' 



h'v 



P 



/2 



i-^lk - '-r + i-^dy + -y ± 



(16) 



(17) 



(18) 



where the positive (minus) sign corresponds to the solution for ip (tp). Transforming to the complex 
coordinates: 

z = X — iy, z = x + iy, (19) 



gives the equation (18) as follows: 



E' 



P 



-Adzdz + jzz + zdz - zdz±l 



At this point, we can define the ladder operators a' and d as: 

which from them, one can write the solution for (p in the following form: 



E' 



2h^v 



2„2 r 



a^a + 1 



(20) 



(21) 



(22) 



where N = a'a, with the eigenvalues n=0,l,..., is the number operator. While for ip, we obtain the 
following solution: 



E- 



,2 2^24 



Here, the two solutions could be packed in one equation as: 



2„,2 



E' 



n^v 



p 



2(n+^)±l 



(23) 



(24) 



At this point one can express the square of the Hamiltonian (4) in terms of the number operator N: 



2 ^ 2h\l f N + 1 
P \ N 

with the following eigenstates and eigenvalues for H: 



^ 



n,K-* 



1 

72 



\n- 1) 
±\n) 



E„ 



fiVp I — 



(25) 



(26) 



It is clear that for n 7^ we have one pairs of eigenstates and eigenvalues but for n = (corresponding 
to E' = 0) we have: 





^ 



0,K+ 



|0) 



(27) 



It is clear that the solution corresponding to the K^ (K^) point^ shows that the probability for the 

electrons to be on the sublattice A (B) is zero. Thus, the zero solution corresponding to K^ (K^) 

implies the localization of Dirac fermions on the B (A) sublattice. In the original experimental paper 

it is argued that since in all the n 7^ energy levels both pseudospin states are filled, whereas in the 

n = level only one is, the density of states in the latter case is 1/2 that of the other levels and 

2 
therefore it contributes only |^ per spin/valley. 

The above argument does not seem entirely satisfactory, since, as it is clear from general solutions 



^The corresponding eigenstate for K is "^o k- 



|0> 




(26), any given eigenstate is normalized to one irrespective of whether one pseudospin component is 

2 

zero or not. Hence the half contribution (f^ ) of the zero energy mode corresponding to the single 
valley index could not be explained in this way because, as we see from the normalized eigenstates, 
electrons localize on just one sublattice, instead of being contributed half on the sublattice A and half 
on the sublattice B. 

There are another explanation for observation of half-integer quantum Hall effect that says since for 
n = we have only one solution (as there is no difference between +|0) and — 10)) the degeneracy of 
this level is half of the other energy levels for which there exist two solutions. What is wrong about 
this conclusion is that existence of just one solution for n = level (and two for others), does not 
simply mean that its degeneracy is twice smaller because one of the two solutions corresponds to the 
negative energy states (holes) and another to the positive energy states (electrons). Therefore this 
assumption could be disregarded. 

Another explanation might be based on this assumption that the level n = is equally shared by 
electrons and holes, meaning that it is half filled with electrons and half with the holes, since there 
is no difference between +|0) and — 10) in this level [10]. In the other words the ground state energy 
level is completely filled with the same types of fermions except the fact that they only differ by their 
charge which does not prevent them from being subject to the Pauli exclusion principle. It is by now 
that one can say the n = energy level contributes l/2(e^//i). Hence, as for the other levels two kind 
of fermions (holes and electrons) with the same number of states contribute in the conductance, they 
contribute twice of l/2(e^//i) per spin/valley. 

The interesting feature that zero energy solutions exhibit, motivate us to seek zero energy modes by 
examining the effect of other types of magnetic fields on graphene's energy spectrum. In fact, when 
the strength of the magnetic field is high the ground state energy level is occupied by more and more 
electrons because the degeneracy of the levels increase and therefore the lowest levels play significant 
role in this case. In this paper, we examine the effect of an inhomogeneous magnetic field which varies 
as -B = (0,0, 1/x) on the lowest energy level and show that it exhibits zero energy solutions which 
is different from those obtained for the case of the constant magnetic field discussed above. As it 
is well-known, this magnetic field occurs around a straight long current-carrying wire (see Fig. 1). 
We first, in the next section briefly discuss the supersymmetric quantum mechanics and the shape 
invariant method [11-14] which turns out to be useful for our investigation. 



2 Supersymmetric quantum mechanics 

One of the methods for solving the quantum mechanical problems is based on finding the relation 
between ground state wave function and the corresponding potential. Considering the Hamiltonian 
H{x) as: 

Hix) = -^ + Vix), (28) 

with associated eigenfunctions and eigenvalues ipnix) and En, respectively, we can write: 

H{x)Mx) = EnM^) (29) 

Now, if one defines Hi{x) as: 



Hi{x)Mx) = H{x) - Eq, (30) 



so that its ground state energy become zero {Eq is the ground state energy of H{x)), we can write: 

It is clear that the two Hamiltonians H{x) and Hi{x) have the same eigenfunctions. Denoting the 
eigenfunctions and eigenvalues of Hi{x) with V'n(^) ^^d .E^, respectively, we can write: 

Hi{x)Mx) = {En - Eo)Mx) -^ Mx)=i^i{x), El = En-Eo (32) 

Now with defining the ladder operators A and A' as: 

A = ^ + Wix), A^ = -^ + Wix), (33) 

ox ax 

where W{x) is called superpotential, one can write Hi{x) in terms of the above operators as: 

Hi{x) = AU = -^ + Vi (x) , (34) 

ax^ 

Here we see that from the relations (31) and (33) one can arrive at the following relation for Vi{x): 

Vi{x) = W\x)-^^^ (35) 

ax 

and keeping in mind that the ground state energy of Hi{x) is zero, we arrive at: 

H, (x) Vo' (x) = i^i^o' (^) = 0, (36) 

which means that A annihilates the ground state wave function ip^{x), i.e.: 

Ml{x) = (37) 

Now it is obvious that from equations 34-36 one can write the ground state wave function with respect 
to the superpotential W{x) and vise versa: 

,o(.).«e- /-•"«- « ,.-,.) ^ 4. n*o(.)^-^^ (38) 

It is by now that we can define Hamiltonian H2{x), partner of Hi(x), which we denote them with 
H^{x) and i7_(x) from now on, respectively, as follows: 

H^{x) = AA^ = -^ + V+{x), (39) 



with 



^2/„, .. ^ , dW{x,ao) 



V+{x) = W\x, ao) + ^^ (40) 

The supersymmetric partner potential V- and V-^ are supposed to be shape invariant if they satisfy 
the following equation: 

V+{x,ao) = V-{x,ai) + R{ao), (41) 



which means that two supersymmetric partner potentials have the same form, but are characterized 
by the different values of parameters a^ and ai . To be more specific, the parameter oi is a function of 
oo, namely, ai = R(ao) with R an independent function of x. Now one can obtain the energy spectrum 
associated to V-{x,ao) simply from the shape invariance condition as: 

n— 1 

E-{ao) = Y,Ria^)^ (42) 

1=0 

where for n > 0: 

E^{ao) = (43) 

while the corresponding wave functions are given by: 

ip~{x, ao) ~ A\x, ao)A\x, ai)...l^(x, ao){x, ai)ipQ{x, an) 

= A\x,ao)i^n-iix^(^i)- (44) 

Note that there are only a few problems that satisfy the shape invariant condition (41). As we show 
in the next section, although the ground state energy level can be obtained analytically, the shape 
invariant condition is not satisfied. 



Zero energy modes corresponding to effect of a varying magnetic 
field 



As we'll show in what follows, graphene spectrum under the influence of a magnetic field which varies 
as inverse of distance, i.e. B = (0, 0, 1/x) exhibits zero energy modes. This magnetic field occurs 
often, as it is the magnetic filed around a long, straight current-carrying wire. In fact because of the 
symmetry of the wire the magnetic lines are circles concentric with it and lie in the planes perpendicular 
to the wire. The magnetic field B is constant on any circle of radius R and is given by: 

where / is the current of the wire and ^o is the magnetic constant. Now if we consider a graphene 
sheet which lies parallel to the axis of wire so that the lines of the magnetic field intersect the graphene 
sheet which is assumed to be in xy-plane, the corresponding vector potential can be written as: 

A={0,qlnx,0), ^B = (0,0,(7-), (46) 

X 

where we have used the Landau gauge and defined q to be: 

At this point, if we go through the same procedure as the case of the constant magnetic field (see 
section 1), in this case, we'll obtain for the S and D (with taking w^? = 1 in our evaluations): 

S = —dx + i{dy — igelnx), 

D = dx + i{dy-iqe\nx). (48) 



In the next step by taking the units such that ?i = c = 1, we can make the foUowing ansatz: 



^(r,t) 



I 1 e'''' 



which leads us to the following equation: 

E^i;{r) = i-V^ + q'^e'^ln'^x + 2iqe\nxdy + —\'il^{i 
We then make the following ansatz for i{j{r) as: 



^(r) 



Jkyy 



fix) 



which plugging it into (50) gives: 



-r^ + k + q e In x — 2qeky In x H 

dx^ ^ X 



fix) = E'fix) 



d^ d^ 

+ Vi{x) = -— + V{x)-eo, 



The above equation is an eigenvalue equation that can be written as: 

H{x)f{x) = efix), e = E^ 
It is by now that we can write the superpotential W{x) in the form: 

Wix) = — alnx H — . 
a 

Then we define Hamiltonian Hi{x) as: 

dx'^ ' ^^ ' dx"^ 

where eo is the ground state energy of H[x) and Vi{x) by the use of the relation: 

V,ix) = W\x)-'^, 

is given by: 

a 6^ 

Vi(x) = a^ln^x-261nxH h ^. 

X a^ 

Now by comparing the two Hamiltonians H{x) and Hi(x), one can get a and b as: 

a = qe, b = qeky, 
which reveals that they are just the same and, therefore, we obtain: 

eo = 0. 



(49) 

(50) 
(51) 

(52) 
(53) 
(54) 

(55) 

(56) 
(57) 
(58) 
(59) 



meaning that the ground state energy level, £"0, is zero. Here we should note that the other energy 
levels can not be derived analytically because the shape invariant condition (41) is not satisfied. We 
also note that for <p{r) the same result is obtained, since we the commutation relation: 



[S,D] = 0, 



(60) 



is satisfied, meaning that the x and y components of dynamical momentum do not commute with each 
other. Thus, we have arrived at a very important result. For graphene under a varying magnetic field 
discussed above, there exists two zero energy modes for which the probability for Dirac fermions to be 
on the sublattice A is the same as sublattice B. However this does not mean that these zero energy 
are different to those obtained for uniform magnetic field in the sense of living electrons on the two 
sublattices. It is because, as we pointed out in the first section, the wave function for the other Dirac 
point are swapped for a constant magnetic field and therefore in an uniform magnetic field electrons 
are present on both triangular sublattices as well as the electrons in a varying magnetic field. The 
whole argue is about the single valley Hamiltonian which One can also obtained foix) as follows: 

/o(x) = N exp(eqxlnx — eqx — kyx) (61) 

where N is the normalization constant and subscribe shows that fo{x) is the wave function asso- 
ciated to the lowest energy level. Note that, here x takes its positive values {x > 0). It may be 
assumed that the existence of the two Dirac points has completely fixed the zero of the energy at 
these points, however, interestingly, for a magnetic field which varies inversely as square of distance, 
i.e. B=(0, 0, A/x^) (where A is a constant), the energy spectrum can be obtained analytically using 
the shape invariant method as [15]: 

I A^e^P 

En = ±VF\ kl ^ (62) 

V (n + l/2 + v/l/4 + Ae(Ae-l))2' ^ ' 

which shows that for ky ^ 0, the n = energy level (-Eq) is not zero, unlike the ground state energy due 
to the effect of the magnetic field from a long current-carrying wire which revealed (even for nonzero 
values for ky ) to be zero. 



4 Implications for experiment 

No zero energy modes are observed when a magnetic field is applied to a system consistent of electrons 
confined in a conventional two dimensional structure, since charge carriers in conventional 2D systems 
obey the schrodinger equation of motion and therefore no massless carriers are imagined for them. 
However, it does not mean that strong magnetic field can not be applied to these systems which give 
rise to the observation of conventional quantum Hall effect. The problem is about the magnetic field 
discussed in the previous section. In fact, one arrives at no analytical solution when the effect of the 
magnetic field B = (0, 0, q/x) (see Eq. (46)) is examined on the massive carries no mater whether they 
behave relativistically or not. This may be the reason why no investigations has been reported up to 
now concerning the effect that this kind of magnetic field might have on conventional 2D system. 
From the above discution we see that graphene could be considered as the only 2D structure that 
investigation regarding the effect of the magnetic field (46) - both from the theoretical and experimental 
point of view - is worth noting. As it is shown in Fig.l, the magnetic lines lie in the planes perpendicular 
to the wire and intersect the graphene's plane. The magnetic field B which is constant on any circle 
of radius R, decrease inversely as the distance increases in the x-direction. 

Here, there is no need to say that the result reported in this paper regarding the existence of zero 
energy modes could be put to the test in contrast with other types of nonuniform magnetic fields such 
as that varies inversely as square of the distance x (see equation 62) and those investigated in [16]. 





Figure 1: Left: The lines of the magnetic field from a current carrying wire which in the xy-plane 
decreases as 1/x, lie in the planes perpendicular to the wire and intersect the graphene's sheet. Right: 
The density distribution associated with the fo{x) for currents /i (green line) and I2 (red line) for 
which I2 = 2/2 • As it is expected the width of the wave-packet grows as the current I -and therefore 
the magnetic field- decreases. 
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In the end of this section, we should note that the magnetic field created by a toroid for a special case 
varies inversely as distance as well. However, it is often used to create an almost uniform magnetic 
field in some enclosed area. One can use Amperes law to obtain the magnetic field inside of a toroid 
with N turns of wire as: 

B = ^ (63) 

27rr ^ ' 

where r which is measured from the center of the toroid is the radius of a circle to which the direction 

of the magnetic field is tangent. In fact the magnetic filed is approximately uniform inside the torus, 

if the radius of toroid, r, is very large compared with the cross-sectional radius of it. But for small 

values of r the magnetic field falls off inversely as r. So our results can also hold for this case as well. 



5 Conclusion 



In this work, we examined the effect of a magnetic field varying inversely as distance on the ground 
state energy level of graphene. One important reason for studying this type of magnetic field- apart 
from this fact that it occurs often-is that it is . In fact it is the magnetic field of a long carrying- 
current wire and, therefore, it can be important when it comes to applications of carbon nanotubes 
and nanowires with graphene. We also showed that graphene under the influence of such a magnetic 
field exhibits zero energy modes which is kind of different from the zero energy modes corresponding to 
the uniform magnetic filed (counted for the observation of the unconventional quantum Hall effect in 
graphene). In fact, contrary to the former case, the zero energy solutions associated to the magnetic 
field B = (0,0, g-) do not show the localization of Dirac fermions on just one sublattice but they 
imply that the probability to find electrons on one sublattice, say A, is the same as other one, say B. 
We also discussed the original interpretation of observation of the half-integer quantum Hall effect in 
graphene which does not seem to be complectly satisfactory because the localization of electrons on 
one sublattice does not imply that the density of states due to the n = Landau level is half of the 
others. We also discussed about how the effect of the two kind of magnetic field which varied as 1/x'^ 
and 1/x on the graphene spectrum could lead to the different results. 

In this work we investigated the effect of a the latter case on the massless Dirac fermions of undoped 
graphene, leading to observation of two zero energy modes which, as we said, are different in the sense 
of living the electrons on the different sublattices (per valley/spin). As we pointed out, considering 
the massive relativistic particles no analytical solution for the lowest energy level is obtained and it 
might be the reason that the potential (45) have not been considered up to now. 
In the end, we should note that at the first sight it might seem strange that the localization of charge 
carriers differs for the varying and constant magnetic field. However, by considering the two Dirac 
points we see that one indicates the localization of electrons on B and another on the A sublattice 
and therefore the equivalency of carbon atoms is not broken. 

Another point which is worth noting here is that the magnetic energy levels obtained from the tight- 
binding model agrees well with that calculated from the kp model [17]. 
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